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ABSTRACT 

In this series of papers, we propose a theory to explain the formation and properties of rings 
and spirals in barred galaxies. The building blocks of these structures are orbits guided by the 
manifolds emanating from the unstable Lagrangian points located near the ends of the bar. In 
this paper we focus on a comparison of the morphology of observed and of theoretical spirals 
and rings and we also give some predictions for further comparisons. Our theory can account 
for spirals as well as both inner and outer rings. The model outer rings have the observed R\, 
R[, R2, R'2 and R1R2 morphologies, including the dimples near the direction of the bar major 
axis. We explain why the vast majority of spirals in barred galaxies are two armed and trailing, 
and discuss what it would take for higher multiplicity arms to form. We show that the shapes of 
observed and theoretical spirals agree and we predict that stronger non-axisymmetric forcings 
at and somewhat beyond corotation will drive more open spirals. We compare the ratio of ring 
diameters in theory and in observations and predict that more elliptical rings will correspond 
to stronger forcings. We find that the model potential may influence strongly the numerical 
values of these ratios. 
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1 INTRODUCTION 



The discs of barred galaxies very often have spectacular sub- 
structures such as grand design spirals, or inner and outer rings. In 
fact, these sub- structures are observed much more often in barred 
than in non-barred disc galaxies. Their formation is intimately 
linked to the evolution of disc galaxies and also gives important 
information on the corresponding underlying potential and the bar 
pattern speed. Therefore many efforts, both the oretical and numer- 
ical, have been made to explain them (see e.g. Linl 1967 . iToomre 
1977L lToomrelll98lL lAthanassoulal[l984 lAthanassoula & Bosma 



1985 and 



Buta & Combesll 19961 for reviews). None of these theo 



ries, however, could explain all observed properties and, particu- 
larly, none could explain both spirals and rings. 

For this reason, we propose in this and in previous papers 
a new theory to explain all these features in a single theoret- 
ical framework. It relies on the dynamics of the unstable La- 
grangian points located near the ends of the bar and of the as- 
sociated manifolds. These guide the motion of the chaotic or- 
bits that constitute the building blocks of the spirals and rings. 
Our theory an d its applications are descri bed in five papers. 
The first two jRomero-Gomez et al \ [2006, hereafter Paper I; 
iRomero-Gomez et al. 1 120071 hereafter Paper II) contain most of 
the theoretical groundwork. They describe in detail the dynamics 
of the region around the two unstable Lagrangian points and of 
the corresponding manifolds. A theory, however, can be dynami- 



cally correct but still irrelevant to a particular application. Thus, 
the aim of the remaining three papers is to check whether our 
theory is applicable to observed spirals and rings. In Paper III 
(lAthanassoula, Romero-Gomez & Masdemontl 12009k we consid- 
ered the orbits guided by the manifolds as potential building blocks 
for spiral and rings and derived their relevant properties. In this 
paper, the fourth of this series, we compare the morphology of ob- 
served and of theoretical spirals and rings to test whether there is 
agreement. We also make a number of predictions based on our 
theory, which can be tested with future observations. In Paper V 
(Athanassoula et al 2009b, to be submitted), the last of this series, 
we will present further comparisons with observations. These last 
three papers can be considered as one entity and therefore have the 
same main title ("Rings and spirals in barred galaxies"), but we will 
always consider them as part of the whole series and refer to them 
as Papers III, IV and V, respectively. 

The comparisons between our theory and the observations 
have, by necessity, to be qualitative. Indeed, in all our calculations 
we have used three specific families of analytic potentials. These 
can not cover all possible forms of observed barred galaxy poten- 
tials, nor are we sure that they share all their relevant properties. 
Furthermore, we do not take into account self-gravity, i.e. we ne- 
glect the potential of the structure itself. If a given quantity de- 
pends considerably on the properties of the potential, we may not 
find good quantitative agreement between our theoretical results 
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and observations, even though our theory is correct and applicable. 
For example, if the ring extent depends strongly on the potential 
used, then one can expect discrepancies between the observed and 
the theoretical extent values, if the adopted potential does not match 
observations in the relevant regions. 

Section [2] contains a brief reminder of the main aspects of 
our theory. This is necessarily simplified as well as condensed, 
so we refer to an interested reader to Papers I, II and III and to 
Romero-Gomez et al. 120091 for more information, a more thorough 
description of the theory and references to other related work. Nev- 
ertheless, the main theoretical concepts and definitions are given 
so as to make this paper as stand-alone as possible, while avoiding 
repetition. In particular we introduce here the Lagrangian points, 
the invariant manifolds and the properties of the associated orbits. 
Section [3] compares the morphological types of our theoretically 
predicted density enhancements with the observed ones. Section |4l 
focuses on the properties of the theoretical and the observed spi- 
rals, such as their sense of winding, their number of arms and their 
pitch angle. Section[5]is devoted to ring properties, such as the fre- 
quency of the various types of rings, the axial ratios of inner and of 
outer rings and the ratio of outer to inner ring major axes. Further 
comparisons with observations, including pattern speed prediction, 
the bar shape, kinematics, radial drift and abundance gradients will 
be discussed in Paper V, the last of this series. There we also give a 
global discussion of all the results in this series. For this reason we 
summarise only briefly the results of this paper in Sect. [6] 



2 THEORETICAL REMINDERS 

Our theory is largely based on the dynamics of the Lagrangian 
points L± and L2 of a two-dimensional barred galaxy system. 
These are located along the direction of the ba r major axis and are, 
in the standard case, unstable saddle points feinnev & Tremainel 
2008). Each of them is surrounded by a family of periodic orbits, 
called Lyapunov orbits dLvapunovll 1 949h . Since these are unstable 
they can not trap around them quasi-periodic orbits of the same 
energ)Q, so that any orbit in their immediate vicinity (in phase 
space) will have to escape the neighbourhood of the correspond- 
ing Lagrangian point. Not all departure directions are, however, 
possible. The direction in which the orbit escapes is set by what 
we call the invariant manifolds. These can be thought of as tubes 
that gu ide the motion of particles of the sam e energy as the man- 
ifolds jGomez et al J 12004 iKoon et al ll200Qh . Each manifold has 
four branches emanating from the corresponding Lyapunov orbit, 
two of them inside corotation (inner branches) and two of them out- 
side it (outer branches). Along two of these branches (one inner and 
one outer) the mean motion is towards the region of the Lagrangian 
point (stable manifolds), while along the other two it is away from 
it (unstable manifolds). We need to stress that the terms 'stable' 
and 'unstable' do not mean that the orbits that follow them are sta- 
ble and unstable, respectively. In fact all the orbits that follow the 
manifold are chaotic, but they are in a loose way 'confined' by the 
manifolds, so that they stay together, at least for a few bar rotations, 
in what could be called a bundle. We propose that these manifolds 
and orbits are the building blocks of the spirals and rings. 



1 We will all through this paper loosely call 'energy' the numerical value 
of the Hamiltonian in a frame of reference co-rotating with the bar, i.e. a 
frame in which the bar is at rest, and denote it by E j . 



These manifolds do not exist for all values of the en- 
ergy, but only for energies for which the corresponding Lya- 
punov periodic orbit is unstable. This means energies within a 
range starting from the energy of the L± or L2 and 
extending over a region whose ext ent depends on the model 
dSkokos. Patsis & Athanassoulal[2002b . 

In Papers II and III we classified the manifolds and the corre- 
sponding orbits in three categories, the homoclinic, the heteroclinic 
and the escaping ones. Homoclinic manifolds and orbits start from 
the vicinity of either the Li, or the L2 and end near the same La- 
grangian point. An example is given in the left panel of Fig. Q] 
Heteroclinic manifolds and orbits start from either near the Li, or 
the L2 and end near the Lagrangian point located at the other end 
of the bar, i.e. the L2 or Li, respectively. An example is given in 
the middle panel of Fig. Q] Finally, escaping orbits start off near one 
of these Lagrangian points but do not end near either this point or 
the symmetric one at the other side of the bar. Their name comes 
from the fact that they escape the vicinity of the bar and reach parts 
of the galactic disc further out, but this does not imply that they 
escape the system and reach infinity. 

The case described so far has five Lagrangian points in total, 
two unstable and three stable, and is referred to as the standard case. 
The formation of the manifolds, and therefore of the rings and spi- 
rals, depends on the existence of such unstable saddle Lagrangian 
points as described above. We have, therefore, applied our theory 
to barred galaxies. It can. however, be applied also to non-barred 
galaxies with other non-axisymmetric perturbations, provided these 
have unstable saddle Lagrangian points as described above and in 
Papers I, II and III. This can include spiral perturbations and will 
be discussed further in Paper V. 

Another alternative, which is interesting both dynamically and 
morphologically, occurs in barred galaxies when the L\ and L2 are 
stable. This can happen if a sufficiently important amount of mass 
is sufficiently concentrated around the Lagrangian points, as de- 
scribed in Paper III. In such cases, the L\ and L2 become stable, 
while four other Lagrangian points appear along the direction of the 
bar major axis, one on either side of the L\ and L2. These new La- 
grangian points are unstable and thus work like the L\ and L2 in the 
standard case. The intricacies of the manifolds in such cases have 
been studied in Paper III, and we will only mention here that the 
gross global morphology of the system stays roughly unchanged 
for a wide range of parameter values. Other things, however, such 
as the axial ratio of the rings can change, as will be discussed in 
Sect. [5] 

In Papers II and III, we calculated manifolds and orbits in 
three different rigid ba rred galaxy potentials. Model A is ta ken from 
I Athanassoula! dl992ah and has a Ferrers bar dFerrersI 18771) of semi- 
major axis a, axial ratio a/b, quadrupole moment Q m and rotating 
with a pattern speed Q p . The central concentration of this model is 
characterised by its central d ensity p c . Model D has a Dehnen-like 
bar potential toehnenll2000h characterised by a strength parameter 
(e) and a scale length (a). The B W model has a bar w ith a Barbanis 
& Woltjer type potential JBarbanis & Woltiedl 19671) . characterised 
by a strength parameter (e) and a scale length (n). The two latter 
models have the same axisymmetric part as model A. Ferrers bars 
have a realistic density distribution, but their non-axisymmetric 
force beyond corotation decreases strongly with radius. The other 
two potentials are ad-hoc, i.e. are not associated with a particular 
density distribution, but can have strong non-axisymmetric forces 
beyond corotation. As in Papers II and III, we use them to represent 
forcings not only from bars, but also from spirals, oval discs and tri- 
axial haloes. A full description of these potentials can be found in 
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Figure 1. Examples of homoclinic (panel a), heteroclinic (panel b) and escaping (panel c) orbits (black thick lines) in the configuration space. In red lines, we 
plot the unstable invariant manifolds associated to these orbits, while in green we plot the corresponding stable invariant manifolds. In dashed lines, we give 
the outline of the bar and, in panels (b) and (c), we plot the zero velocity curves in dot-dashed lines. 



Appendix A of Paper III. These bar potentials will be used here 
together with different types of rotation curves. As in our previous 
papers, we measure angles in the standard mathematical sense, i.e. 
starting from the positive x axis and increasing anticlockwise. 

As shown in Papers II and III, the properties of the manifolds 
depend strongly on the relative strength of the non-axisymmetric 
forcing at and somewhat beyond corotation. We measure this with 
the help of the quantity Qt,L 1 , which is the value of Q t 

Q t (r) = (0$(r, 0)/dO)max/(rd$ /dr), (1) 

at t\l, the radius of the Lagrangian point L\ or L2, i.e. Qt,L t = 
Q t (r = tl). This is not the same as the strength of the bar, which is 
often measured by Q b, i.e. the maximum of Q t over all radii shorter 
than the bar extent. The radius at which this maximum occurs can 
be small compared to vl. So Qb is not necessarily a good proxy for 
Qt,Li and vice versa. Nevertheless, for the sake of brevity, we will 
often replace in our discussions "non-axisymmetric forcings which 
are relatively strong at and beyond corotation" simply by "strong 
bars", or "strong non-axisymmetric forcings". 



3 MORPHOLOGY 

Our aim in this work is to show that the dynamics of the mani- 
folds emanating from the unstable Lagrangian points can explain 
the formation of spirals as well as of inner and outer rings in barred 
disc galaxies. Thus, in this section, we will make the link between 
the morphology of theoretical and of observed spirals and rings. 
At this stage, it is premature to make a detailed comparison for 
one specific galaxy. Instead, we will make global comparisons and 
make sure that there are no major differences between models and 
observations and no obvious discrepancies which would disqualify 
our theory. 

Any candidate theory should reproduce correctly the observed 
spiral and ring morphologies. Our theory can indeed accomplish 
this, as we will show in the next subsections. The top two rows of 
panels of Fig. [20 show four galaxies that contain archetypes exam- 
ples of the characteristic morphologies that have to be reproduced; 

2 The images of NGC 1365 (first column, first row) and of ESO 507-16 
(third column, first row) are taken from the Digital Sky Survey and those of 
NGC 2665 (second column, first row) and ESO 325-28 (third column, first 
row) were reproduced fro m http://bama.ua.edu/~rbuta/devatlas (see also 
Buta, Corwin & Odewahn 2007). The three schematic plots in the second 
to fourth column and second row are from Buta & Crocker (1991) and are 
reproduced by permission of the American Astronomical Society. 



namely spirals, inner rings and the three types of outer rings, Ri , 
R2 and R1R2. They will be described in detail in the following 
subsections. The upper row of panels gives the images and the sec- 
ond row the schematic plots, bringing out the characteristic mor- 
phologies of the galaxies, which the manifolds have to reproduce 
The third and fourth rows give manifold examples which repro- 
duce qualitatively the main characteristic features well - the former 
viewed from an orientation similar to that of the observed galaxy in 
order to permit a better comparison, the latter viewed face-on. 

Note that we have not used the potential of each specific 
galaxy, nor have we made a search for the best fitting manifold. 
We simply found, amongst the very large number of models for 
which we have calculated manifolds, one which roughly represents 
the desired morphology. 



3.1 Spiral structure 

The left column of Fig. [2 corresponds to NGC 1365, a galaxy with 
a strong bar and a clear, well-developed, grand design spiral struc- 
ture. It has two spiral arms, each emanating from one of the two 
ends of the bar and extending well into the outer disc. The lower 
two panels in this column show a model qualitatively reproducing 
this morphology. It also has two spiral arms, one emanating from 
each end of the bar and extending outwards. The spiral winding also 
appears similar and will be discussed in more detail in Sect. 14.31 
Thus the model matches qualitatively the morphological properties 
of the spiral very well. 

A further point to note is that in the case of escaping trajecto- 
ries - i.e. trajectories that wind from L± (or Li) outwards and do 
not end at the same or the other Lagrangian point - the manifolds 
become wider as the distance from the centre increases (see e.g. 
bottom left panel of Fig. [2]>. Under the assumption that no mass is 
created or destroyed, this implies that the density in this arm com- 
ponent should decrease with increasing distance from the centre. 
However, if other arm components - e.g. a weak, leading compo- 
nent - are present, this will lead to bumps on an otherwise smoothly 
decreasing density profile. 



3.2 Inner rings 

Inner r ings have roughly the same size as the bar dButall 19861) . iButal 
dl995h studied their orientation with respect to the bar major axis 
using his catalogue of southern ring galaxies (CSRG). Comparing 
with model distributions of apparent bar/ring position angles, he 
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Figure 2. Ring and spiral arm morphologies. The upper row shows four galaxies, typical of the spiral, R\, R2 and R1R2 morphologies, namely NGC 1365, 
NGC 2665, ESO 325-28 and ESO 507-16, respectively. The second row shows schematic plots of these galaxies, bringing out the spiral and ring structures. 
The third row gives examples from our manifold calculations which have a similar morphology, projected in roughly the same way as the observed galaxy or 
as the schematic plot. Finally, the fourth row gives the face-on view of these manifolds with the bar along the x axis. Note that, contrary to all other cases, the 
sense of rotation in the lowest left panel is clockwise in order to ease the comparison with NGC 1365. 



found that inner rings are preferentially elongated in the direction 
of the bar major axis. 

In our models, some of the inner manifolds give rise to inner 
rings while others contribute to the bar, particularly its outer parts. 
If the Li and L 2 are not too strongly unstable, the inner manifolds 
retrace the same path after one revolution around the galactic cen- 
tre (Paper III). In this case it is easy to get the outline of the ring. 
A good example can be seen in the lower two rows of the sec- 
ond column of panels in Fig. |2] Other examples can be found in 
Papers I, II and III. In these cases, the inner manifold branches out- 
line a ring which has the same morphology as the observed rings. 
Namely, it has roughly the same extent as the bar and is elongated 
parallel to the bar major axis. In cases with strongly unstable La- 
grangian points, the inner manifold branches retrace the same path 
only after a few revolutions. Thus the shape of the corresponding 
density enhancement should be calculated by the self-consistent su- 
perposition of the inner manifold loci, and this is well beyond the 
scope of this paper. To summarise, the inner manifold branches ei- 
ther reproduce well the shape and location of observed inner rings, 
or populate mainly the outer parts of the bar. 



3.3 Outer rings 

Outer rings are considerably larger than inner ones, with a major 
axis of the order of twice the bar size. Depending on the re lative 
orienta tion of the major axes of the outer ring and of the bar, iButal 
i 19951) distinguished three different types of outer rings. Rings of 
Ri type are elongated perpendicular to the bar major axis. Rings 
of R2 type are elongated along it and R1R2 have a part elongated 
along the bar and a part elongated perpendicular to it. 

NGC 2665, ESO 325-28 and ESO 507-16, shown in Fig.0 
(second to fourth column of the top row), are prime examples of 
these three types of outer rings. Note that NGC 2665 includes a 
characteristic feature of R± rings. Namely at or near the direction of 
the bar major axis, the ring approaches the ends of the bar. The fea- 
tures thus forme d are called dimples feutalll986l : iButa & Crocker! 
ll99ll : lButall995h and give to these outer rings a characteristic shape 
reminiscent of the number 8. Their schematic plots are given in the 
second row. The third and fourth rows show examples of manifolds 
with similar morphological characteristics. 

The example given in the bottom two rows of the second col- 
umn reproduces well the characteristic morphology of rRi rings in 
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general and of NGC 2665 in particular. It has both an inner and an 
outer ring. The latter is elongated perpendicular to the major axis 
and has the characteristic morphology of a figure '8', with clear 
dimples near the two ends of the bar. Thus, the manifolds reproduce 
well all the main observed morphological features of observed rR\ 
rings (see also Paper I). 

The example in the third column has a different morphologi- 
cal type; that of an R 2 ring. There is a spiral starting from the end 
of the bar, which, after winding by roughly 90°, joins with an outer 
ring. The latter is elongated along the direction of the bar major 
axis, but its ellipticity is rather small, i.e. it is not too far from cir- 
cular. Also it does not have the dimples seen in the R± rings. These 
morphological characteristics are shared by both the galaxy and the 
manifold. 

The example in the fourth column is of the Ri R2 morphology 
and has two outer structures, i.e. two outer rings. The inside one is 
elongated perpendicular to the bar major axis and has characteris- 
tics dimples near the ends of the bar, i.e. it is of the Ri variety. The 
outermost one is elongated along the bar major axis, has a relatively 
small elongation and no dimples in the direction of the bar major 
axis. All these morphological characteristics are clearly seen in the 
manifold examples in the last two rows. Also, in all the examples 
of Ri R2 morphology we have calculated so far, the Ri feature lies 
within the R2 feature and this also holds for all the observed cases 
we have seen. 

To summarise, we can say that the main morphological fea- 
tures of the three types of outer rings are qualitatively very well 
reproduced by our theoretical results. 

3.4 Bar strength and morphological types 

In Paper III we showed that the morphological type - i.e. whether 
we have a spiral, an R±, an R2, or an R1R2 ring or pseudo-ring 
- depends on the strength of the non-axisymmetric forcing at the 
corotation radius. Relatively weak non-axisymmetric forcings give 
Ri and R[ type morphologies, while stronger forcings give spirals 
and the remaining types of rings. 

A related result was found in Paper II, where we showed that 
the morphology can also depend on whether the rotation curve is 
rising, flat or slightly falling in the relevant part of the disc, i.e. at 
and somewhat beyond the Lagrangian radius. In examples given in 
Fig. 8 of that paper, we showed that, for the same bar, a rising ro- 
tation curve can give a Ri or R[ morphology, while a decreasing 
rotation will give an R2 or a spiral morphology. This can now be 
understood from the results presented in Paper III. Indeed a differ- 
ent axisymmetric model will change the values of Qt at and beyond 
the Lagrangian radius, and therefore the morphology and the prop- 
erties of the rings and spirals. Rising rotation curves imply more 
mass in the relevant part of the galaxy and therefore a stronger ax- 
isymmetric force and a smaller Qt,L ± value, the opposite being true 
for declining rotation curves. Thus the result found in Paper II is in 
fact a particular aspect of what was found in Paper III, namely that 
the morphological type depends on Qt,L 19 the value of Q t at and 
somewhat beyond the Lagrangian radius. 

This is a clear prediction, which can be tested by observations. 
In doing so, two points need to be kept in mind. The first one is that 
the relative importance of the non-axisymmetric forcing depends 
not only on the bar or spiral, but also on whether the rotation curve 
is flat, rising or somewhat falling (Paper II). Thus such a test is not 
possible solely with photometric data, but necessitates also kine- 
matics. The second point to be kept in mind is that a spiral forcing 
can produce spirals easier than a bar forcing of comparable ampli- 



tude. Thus, when comparing spiral to ring morphologies one should 
keep in mind that the potential of the spiral has to be also taken into 
ac count. This has been already argued, albeit i n a different context, 
bv lLindblad. Lindblad & Athanassoulal dl996h . 

Making a quantitative comparison is, at this stage, not straight- 
forward because our theory is not self-consistent and because we 
only use a bar forcing, while in real galaxies the rings and spirals 
are formed under the combined effect of the bar and of their own 
pote ntial. The bar forci ng is of course much stronger than the spiral 
one. iButa et al\ d2Q05h give for the SB galaxies in their sample for 
the bar an average Qb of 0.294 ± 0. 21 and for the spiral an av- 
erage Q s of 0.171 ± 0.014. Similarly. IBurbala et all (120091) . using 
a sample of isolated barred galaxies of type Sb-Sc, find 0.243 ± 
0.015 and 0.166 ± 0.009 for the mean values of Qb and Qs, re- 
spectively. For the medians they find 0.222 and 0.165, respectively. 
These numbers can not be directly translated into a ratio of relative 
effects, because the maximum of the spiral forcing occurs at larger 
radii than that of the bar forcing. They nevertheless clearly show 
that the effect of self-consistency will change somewhat the limits 
between the zones with different morphologies (R±, R[, spirals, 
R2, R1R2 etc; see Paper III), shifting them to lower bar strengths. 
The amplitude of this shift will depend on the strength of the spiral, 
so t hat no quanti t ative c omparison is possible at this stage. 

Buta et al. ] J2005h measured the bar strength for 147 spi- 
ral galaxies in the statistically wel l-defined Ohio State U niversity 
Bright Galaxy Survey fOSUBGS: lEskridge et al\ fcOOOh ) sample 
and found that the average bar strength in barred galaxies with spi- 
rals (SB(s) types) is < Qb >— 0.329, i.e. higher than the average 
for all barred galaxies, which is < Qb >— 0.294. They also find 
that the average bar strength in barred galaxies with inner rings 
(SB(r) types) is < Qb >= 0.201, i.e. lower than the average for 
all barred galaxies. These trends are rough, but still are in good 
agreement with our theoretical predictions from Paper III. 



4 MAIN SPIRAL PROPERTIES 
4.1 Why are spirals trailing? 

Simply by viewing a spiral galaxy projected on the sky, it is 
not possible to say whether the arms are trailing or leading, 
even if one knows the sense of rotation of the mat e rial w ithin 
the galaxy disc. As explain ed by iBinnev & Tremainel d2008h and 
iBinnev & Merr ifield ( 1998|), one needs also to know which side of 
the galaxy is nearest to us. This can be obtained b y s tudying the ob- 
scuration by dust (see e.g. ldeVaucouleurslll959l and lSandagelll96ll 
for such pioneering studies). In almost all cases where the result is 
unambiguous the arms are found to be trailing. In our own Galaxy 
also the arms are trailing. It is thus generally admitted that spiral 
arms in disc galaxies are trailing. How does that compare with our 
theoretical predictions? 

The anti- spiral theorem jLvnden-Bell & O striked! 1967b states 
that if a steady- state solution of a time-reversible set of equations 
has the form of a trailing spiral, then there must be an identical 
solution in the form of a leading spiral. Thus, density wave theory 
can not distinguish between trailing and leading solutions, unless 
one takes into account the growth of the spiral. 

Swing amplified s pirals, however , are not stationary, so do not 
have this shortcoming jToomrelll98ll) . The inwards going wave is 
leading and will be reflected at the centre of the galaxy into an 
outwards-going trailing wave which will be amplified in the corota- 
tion region. Therefore, according to this theory, both a trailing and 
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a leading component will be present, but since a very strong am- 
plification occurs, the trailing component will have a much larger 
amplitude than the leading component. Thus, visually one will see 
a trailing spiral, whose amplitude will present oscillations with ra- 
dius due to the interaction of the leading and tr ailing spiral phases 
dAthanassoulalll984[l992bl : IPuerari et al.ll2Q0Qh . 

Since the potentials we have used so far are stationary in the 
rotating frame of reference, the anti- spiral theorem will apply. In- 
deed, if one reverses the arrow of time in the case of the unstable 
manifolds - which have a trailing spiral shape - one finds the stable 
manifolds - which are of the same shape but leading. This is true 
whenever the bar does not evolve. Material, however, needs ini- 
tially to get trapped into the manifolds, and this can only be done 
while the bar is being formed, or while it grows. This introduces 
time evolution and thus breaks the degeneracy between trailing and 
leading arms. To understand this further, we made a number of re- 
sponse simulations. For this we considered 10 000 particles initially 
in circular orbits around the galactic centre and grew gradually (in 
one bar rotation) a bar in the disc. 

Let us start by describing the results of a case which, when 
the bar has reached its final amplitude, will have a spiral manifold 
morphology. The response is a trailing spiral which forms grad- 
ually while the bar grows. This means that at any time there are 
more particles in the trailing arm than in the background or the cor- 
responding symmetric leading arm. This, however, does not, on its 
own, prove that the unstable (trailing) manifolds trap particles as 
the bar grows, because particles could still traverse the arm, which 
would then be a density wave. 

To reveal the differences between the unstable and stable man- 
ifold and to prove that it is mainly the unstable manifold that traps 
particles, we need to work in phase space rather than configuration 
space, i.e. we need to consider velocities as well as positions. We 
thus chose a number of simulation particles at random in the radial 
range of the spiral arms and in the energy range where the stable 
and unstable manifolds exist. For each of these, we calculated the 
manifolds at the particle energy and found the tube they outline in 
phase space. We then checked whether the particle is trapped by the 
stable or the unstable manifold, or whether it is not trapped at all. 
We repeated this for several time steps and a large number of parti- 
cles and from these results we can get some statistical information. 
We find that of the order of 30 - 50% of the particles are trapped by 
the unstable manifold and only - 2% by the stable ones. The re- 
maining particles are either not trapped by any manifold, or it was 
not possible to give an unambiguous answer due to their locations 
(e.g. close to the Lagrangian points). This very important numeri- 
cal difference between the stable and the unstable manifolds shows 
that our theory can make a strong prediction concerning the sense 
of winding of the spiral arms. Namely, it predicts that as the bar 
grows in a galaxy, the unstable manifolds will trap mass elements. 
Since the unstable manifolds have a winding in the trailing sense, 
this means that our theory predicts that spiral arms generated by the 
manifolds of the L± and the h Lagrangian points should be trailing, 
as observed. 

The situation is different in the case of a bar which, after the 
bar has reached its maximum amplitude, has manifolds with the 
shape of rings. As we already discussed in Paper II, in this case 
the unstable and stable manifolds connect in phase space, so it is 
not possible to say whether a particle belongs to the stable or the 
unstable manifold. Thus, in this case, both stable and unstable man- 
ifolds can be populated, so that there is no preference of trailing or 
of leading windings. This is in good agreement with the observed 
shape of rings in real galaxies, which have no preferred winding. 



Thus to summarise, in the case of spirals, material gets trapped 
mainly in the manifolds as the bar grows, mainly in the unstable 
one. This accounts for the preponderance of trailing (as opposed to 
leading) spirals. In the case of rings, the stable and unstable mani- 
folds communicate in phase space, so there is no preferred sense of 
winding. 

4.2 Number of arms 

In the vast majority of cases, spirals in barred galaxies have two 
spiral arms. This is in good agreement with our theory which re- 
lates the spiral arms to the unstable Lagrangian points. In the stan- 
dard case, a barred galaxy has two unstable Lagrangian points, L± 
and L2, which are saddle points (Paper I). Each one accounts for 
a spiral arm. Thus our standard theory predicts that barred galax- 
ies will have two-armed grand design spirals emanating from a lo- 
cation near one of the two opposite ends of the bar. Since all the 
forcings we have used so far are bisymmetric, we found by neces- 
sity two symmetric spirals. In real galaxies, however, there are al- 
ways asymmetries in the forcing potential, and this will lead also to 
asymmetries between the two spirals, or to small angular shifts of 
their beginning with respect to the tip of the bar. There is thus very 
good agreement between the results of our theory and observed spi- 
rals regarding the number of spiral arms. 

Although rare, there are, nevertheless, a few observed barred 
galaxies wi th more arms, in particular with a four armed structure 
dButalll995h . A good example is ESO 566-24. Can our theory ac- 
count for such structures? 

If more than two Lagrangian points with saddle behaviour 
were present, our theory could predict spirals of higher multiplic- 
ity. The most straightforward way of achieving this is to include 
a cos(mO) forcing, with m > 2 and a sufficiently high ampli- 
tude. Making a full analysis of such potentials, of their origin and 
of the corresponding manifolds is beyond the scope of this paper. 
We will, nevertheless, examine the iso-effective-potential curves 
and the type of stability of the corresponding Lagrangian points 
for some specific cases, in order to get some insight and to test 
whether such structures are possible. We will also, for simplicity, 
limit ourselves to m = 4 cases, but the extension to other m values 
is straightforward. 

Ferrers bars have eightfold symmetry and forcing terms m = 
2 and higher. Increasing the bar mass increases both the m = 2 and 
the m = 4 terms leaving their ratio unchanged. Increasing the bar el- 
lipticity, however, changes the m = 4 over m = 2 ratio; thinner bars 
having larger val ues of this ratio. In su ch bars the L4 and L5 can 
become unstable jAthanassoulalll992aL but they become complex 
unstable, not saddle points. These do have manifolds, both stable 
and unstable, but it is unclear whether these can account for spiral 
arms and at a first glance it does not seem very promising. Further 
theoretical work is necessary to clarify this issue. 

The simplest and most straightforward bar potential having 
the desired properties is of the form 

$ b (r, 0) = A 2 (r) cos(26>) + A 4 (r) cos(46> + 0(r) + O ), 

This is bar-like within corotation with only an m=2 forcing and 
spiral-like with m = 4 outside it. Fig. [3] shows a specific example, 
whose iso-effective-potential plot looks promising. Around corota- 
tion there are four unstable Lagrangian points of the saddle type, 
similar to those at L\ and L2 in the standard case. At a somewhat 
larger radius there are four stable Lagrangian points, similar to the 
L4 and L5 in the standard case. Introducing such an m forcing is 
not the only way to obtain four unstable Lagrangian points with a 
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Figure 3. Iso-effective contours for a potential with an m = 2 symmetry 
inside corotation and an m = 4 outside it. There are four Lagrangian points 
in the region around corotation and four more somewhat beyond that. The 
former are of saddle type and are marked with a red star, the latter are stable 
and are marked with a red dot. The four-armed barred galaxies, such as ESO 
566-24, could be due to potentials of similar type. 

saddle point behaviour. This can also be achieved with an m = 2 
only, provided it has a non-monotonically decreasing profile of the 
appropriate shape. Again we stress that these are very simple ad hoc 
potentials, which may not necessarily be linked to a realistic den- 
sity distribution and do not represent any galaxy in particular. It is 
nevertheless useful to consider them, because they may give some 
insight as to wh at kind of potential can explain the spiral structure 
in ESO 566-24 dRautiainen. Salo & Buta1l2004l) . or other, undoubt- 
edly very rare, galaxies of this type. 

To summarise, our theory predicts that the vast majority of 
barred spirals will be two armed, but that, for certain types of po- 
tentials, it is possible to have more arms. This is in good agreement 
with what is indeed found by observations. 

4.3 Spiral arm shapes 

We showed above that the manifolds have shapes that can re- 
produce spiral arms. But how does the winding of the mani- 
fold loci compare with that of spiral arms in real barred galax- 
ies? To answer this question we deprojected the image of NGC 
1365 (upper two panels of the left column of Fig. [3, using view- 
ing angle values obtain e d by r ounding off the estimates given by 
Ijorsater & van Moorsell QJ995), namely a position angle of 220° 
and an inclination angle of 50°, and we measured at each radius 
the coordinates of the maxima along the arm. We can now compare 
the loci outlined by these measurements to the loci of a typical spi- 
ral manifold. For this we chose a D model with e = 0.05 and tl 
= 5. In Fig. HI we can see that the comparison is very satisfactory. 
In particular, both observed spirals and calculated manifolds share 
the common property that after reaching a maximum radius they 
fall back towards the inner parts of the galaxy. For a more precise 
comparison we would need to calculate the potential of NGC 1365 
from some image in the near infrared as well as kinematics, but this 
is beyond the scope of this paper. 

Thus our theory, explaining the formation of spiral arms in 
barred galaxies with the help of the invariant manifolds of the L\ 
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Figure 4. Comparison of the manifold loci on the 6 = ln(r) plane (blue 
on the online version, black in the paper version) with the spiral arm loci in 
NGC 1365 (x signs, green in the online version, gray in the paper version). 
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Figure 5. Pitch angle i as a function of bar strength measured by Qt,L 1 - 
We plot the results only for three values of r\ - namely tl = 5, tl = 5.5 
and tl = 6 - but we have verified that other realistic values also show the 
same trend. We include in the plot models from all three families, A, D and 
BW. For A models we include in the plot only models with a Ferrers' index 
n = 1, but we verified that models with n = lie on the same curve. 

and L2 Lagrangian points, explains the characteristic winding of 
the spiral arms observed in a large fraction of barred galaxies. 
Namely, as the angle increases the radius first increases and, upon 
reaching a maximum radius, 'falls back' towards the bar region, 
or towards the region of the other spiral. This property has not 
been so far reproduced with density wave models, where the spi- 
ral arms tight en up as they reach out er Lindblad resonance, but do 
not fall back (ILin. Yuan & Shull 19691 and later works of this group, 
lAthanassoulalll98ol for an approach not relying on the WKBJ ap- 
proximation). 

4.4 Trends with the pitch angle 

In order to assess how various bar parameters affect the winding 
of the spiral, we obtained the value of the pitch angle by fitting a 
logarithmic spiral 

— 60 — cot i ln(r/ro) 

over the main relevant part of the manifold, for many models. We 
examined each of our three model families separately (families A, 
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D and BW), and used as a measure of the bar strength Q m , e and 
e, respectively. We find that (not shown here), for all three fami- 
lies, the pitch angle increases with increasing bar strength and the 
form of the dependence is similar. In other words, the spirals cor- 
responding to stronger bars are more open than those from weaker 
bars. This is a direct consequence of the spiral properties described 
in section 3 of Paper III. To our knowledge, no such relation has as 
yet been reported by observations. It would be thus very useful to 
test whether such a relationship can indeed be found. In order for 
this to be possible, we need first to check whether the trend persists 
when the pitch angle of the various models is plotted as a function 
of a common bar strength measure. For the latter we use Qt,L 1 and 
show the results in Fig. [5] for the three model families and three 
pattern speed values. The trend is clearly present, with relatively 
little spread due to the different models. The pitch angle increases 
steeply as Qt,L 1 increases from to 0.2, and has a more shallow 
increase beyond that. It is not easy to confront this result with mea- 
sure ments from the barred galaxies of the O SU and the Frei sam- 
ples feskridge et al Il200d : iFrei et al. Ill996h for mainly two rea- 
sons. The first difficulty lies in the deprojection of the galaxies, 
since the value of the position and inclination angles can be heavily 
influenced by the outer par ts of the arms. This is the well known 
Stock's effect dStocklll955h and may influence tightly wound spi- 
rals differently from open ones. This could introduce a bias in the 
analysis and thus has to be treated with particular care. The second 
difficulty comes from the fact that the quantity Qt,L 1 is not easy to 
obtain from observations since it implies knowledge of the corota- 
tion radius and some information of the dark matter content in this 
region. 

Such a relation has not been found by density wave theory be- 
cause this has been fully developed only in the linear regime, so that 
a change in the bar am plitude will on ly c hange the amplitude of the 
spiral (for reviews, see lToomrell 19771 and lAtha nassoula 1984). It is, 
however, present in the response simulations of lSchwarzl dl984h . as 
expected, since, as we showed in Paper III the gas-response spirals 
and rings obey manifold dynamics. 



5 MAIN RING PROPERTIES 
5.1 Frequency of ring types 

iButal dl995h classified 3692 galaxies in his CSRG catalogue and 
found that about 53% of them have an outer ring. 25% of the galax- 
ies with an outer ring have an Ri ring or R[ pseudo-ring morphol- 
ogy. About 15% have an R2 ring or R' 2 pseudo-ring and 4% an 
Ri R2 ring morphology. The rest of the ringed galaxies have a sin- 
gle near circular outer ring. 

Our model reveals what part of the parameter space corre- 
sponds to what type of rings (see Paper III), but we can not translate 
this into frequencies of observed ring types because our models are 
not exact replicas of observed bars. It is possible, however, to make 
a few relevant comments. In particular, in Fig. 4 of Paper III we ex- 
amined the morphology of models with Ferrers bars, which are, to 
date, the most realistic analytic bar model potential. For these, we 
found that in most cases the manifold morphology is of Ri , or R[ 
type. This is in good agreement with the fact that Ri is the most fre- 
quently observed of the three ring types. Although this agreement is 
very welcome, it should not be considered as a confirmation of our 
theory. On the other hand, it would have been worrying if observa- 
tions showed that R2, or R1R2 rings were most common since this 
would have implied a rather bizarre distribution of bars on the (bar 



strength, bar pattern speed) plane, and/or very strong differences 
between Ferrers potentials and the potentials of observed bars. 



5.2 Ratio of ring sizes 

5.2.1 General remarks 

As shown in Paper III, ring sizes depend strongly on the model po- 
tential used. Moreover, our models are not exact representations of 
observed bars and we neglect the self-gravity of the ring structure. 
Therefore, quantitative comparisons of model and observed ring 
ratios may well not show agreement. 

There is a further complication concerning the measurement 
of inner ring sizes in models. As described in Paper III, weak bars 
have weakly unstable L\ and L2 Lagrangian points and, after a 
revolution around the galactic centre, the inner manifold branches 
retrace the same path, thus forming the inner ring. In these cases it 
is relatively easy to measure the length of the major and minor ring 
diameters. On the other hand, strong bars have strongly unstable L\ 
and L2 and the inner manifold branches cover a new path after the 
first revolution around the galactic centre, retracing the same loci 
only after a couple or a few revolutions (Papers III and V). One of 
these loci could correspond to the ring, and the others to the outer 
parts of the bar, or all could correspond to the bar. In the latter case 
there will be no inner ring, while in the former one it is not easy to 
find the ring outline. In principle, this should be measured from the 
self-consistent superposition of the corresponding manifold parts. 
This, however, is beyond the scope of this paper, so in what fol- 
lows we measured, whenever relevant, the size and shape of the 
inner ring from the loci of the first revolution of the inner mani- 
fold branches. This is arbitrary and thus could lead to considerable 
systematic errors. In particular, it could lead to a considerable un- 
derestimate of the ring minor axis. This will be discussed further 
below. 

Further reasons for discrepancies can stem from the observa- 
tions. Galaxies are observed as projected on the sky, i.e. are not 
necessarily face-on, and thus the axial ratios can not be directly 
compared to models. IButal (1 19951) gave statistical estimates of the 
deprojected axial ratios from a careful analysis using a model and 
relying on assumptions. The obvious and indisputable assumption 
is that the observed galaxies are oriented randomly with respect to 
the line of sight. Nevertheless, since rings can not be seen in edge- 
on, or near edge-on, discs, there is a bias against highly inclined 
cases, which is dealt with in the modelling by an approximation. 
Further assumptions include the fact that both discs and bars are 
vertically infinitesimally thin, that the ring shapes are well approx- 
imated by ellipses, and that the intrinsic shape of the axial r atio 
distributio n is either t op-hat with a narrow extent dButalll986K or 
Gaussian feutalll995l) . Moreover, problems may also come from 
selection effects in the catalogue, especially for rings with small 
diameters. Thus, the errors stemming from observations can also 
be important. They will, moreover, be more severe for inner than 
for outer rings and will b ring an arti ficial reducing of the number 
of small axial ratio cases teutalll995l) . 

Despite the above arguments, we attempt in this section not 
only qualitative but also quantitative comparisons and estimates, 
but it should be kept in mind that the latter may be affected by 
non-negligible error, particularly in the case of inner rings. Let us 
call D and d the lengths of the major and minor axes, respec- 
tively, of the outer ring and by Di and di the same quantities for 
the inner ring. The relevant quantities which we can use to make 
comparisons with observations are d /D , di/Di and D Q /Di. 



Rings and spirals in barred galaxies. II. Ring and spiral morphology. 9 




-6 D 6-6D 6-6D 6 -6 D 6 —5 D E —5 D E -6 D t -6 D 6 — O D fl 



r L =E>, r L =6,35 r L =5,5 r L =5,75 r t =6, r L =6,25 r L =6,6 r L =6,75 r t =7. 

Figure 6. Effect of the bar strength and of the Lagrangian radius r l on the manifold loci for a model with a Ferrers bar and a slightly decreasing rotation 
curve (see text). The strength of the bar is parametrised by the quadrupole moment Q m - Each panel corresponds to a pair of (Qm, values, given on the 
left (in units of 10 4 ) and on the bottom of the figure, respectively. The manifolds are plotted in a colour which is determined by their morphology: green for 
R±, orange for R' x and black for spirals. The thin, black lines give the outline of the bar. Compared to fig. 4 of Paper III, which has an identical bar model, this 
plot shows how a change in the axisymmetric component can affect the morphology of the response and the axis ratios of the rings. 



Table 1. Results for models with lenses: Axial ratio of the inner and outer 
ring sizes (di/Di and d /D , respectively) and Qt,L x as a function of the 
relative lens mass, fi ens (see text for more information). 
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Table 2. Results for models with unstable or stable L\ and L2 Lagrangian 
points: Axial ratio of the inner and outer ring and stability of the L\ and L2 
as a function of the relative stabilising mass F (see text and Paper III for 
more information on the models). 
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5.2.2 Further models 

Since the ring shapes and the ratio of ring major axes diameters 
may depend crucially on the adopted potential, we will, in the next 
few subsections, consider not only models A, D and BW (see Paper 
III), but also other potentials. 

To model galaxies with lenses, we add an extra Ferrers com- 
po nent which corota tes with the bar and has the same major axis as 
it dKormendvl 1 1 979h . The remaining parameters were taken as for 
the fiducial model of Paper III, i.e. a — 5, a/b — 2.5, n = 6, 
Q m — 4.5 x 10 4 , p c = 2.4 x 10 4 and n = 1, all measured in 
the units given in Appendix A of that paper. The lens has an axial 



ratio b/a = 0.9 and a mass Mi ens = fiensMb, where Mb is the 
mass of the bar. The results are given TableQ]for two such series of 
models. In the first one the total mass of the model was kept con- 
stant (columns two to four, marked 1), so that a certain fraction of 
the bar mass was converted into a lens. For the second one the lens 
was simply added on to the model, without any change of the bar 
mass (columns five to seven, marked 2). For the first series, increas- 
ing the lens fraction decreases Qt,L 1 (as expected), and renders the 
inner and outer rings more axisymmetric. For the second series of 
models, adding a lens makes a less clear-cut effect, because the to- 
tal mass of the model does not stay constant. It is, however, clear 
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that it makes the inner ring more axisymmetric and does not have a 
big effect on the shape of the outer ring. 

We also considered models in which the L\ and L2 La- 
grangian points are stable. As in Paper III, this stability was 
achieved by adding t wo mass concentrations in th e form of 
Kuz'min/Toomre discs jKuz'mirJ[l956l : lToomrd ll963) each cen- 
tred on one of the L\ and L2 Lagrangian points. In the examples 
given here we adopted for these discs a scale-length r a — 0.6 
and we tried several values of the fractional mass F, where F = 
Mkt(t < r a )/Mb(r < r a ) is the ratio of the mass of one of the 
Kuz'min/Toomre disc (Mkt) to that of the bar (Mb), both calcu- 
lated within a radius equal to r a from their respective centres (see 
Sect. 5 in Paper III for more information). These values, the result- 
ing stability and the measured outer ring axial ratios are listed in 
Table [2 We note that, as F increases past a certain threshold, the 
Li and L2 turn from unstable to stable. An increase in F makes the 
inner ring considerably more elongated, and has a smaller effect on 
the outer ring axial ratios, which furthermore is not always mono- 
tonic. As a result, the ratio of outer to inner major axes increases 
considerably with F. 

We also analysed the results for two further models, similar 
to model A but with a different axisymmetric component. For the 
first one we chose the axisymmetric potential leading to the circular 
velocity curve denoted by D in Paper II (see Fig. 7 in that paper.). 
This is slightly decreasing in the regions around corotation and im- 
mediately beyond it. Thus, Qt,L ± will be somewhat larger than for 
model A. We used the same bar model and calculated the manifolds 
for the same grid of parameter values as for model A in Paper III. 
We call this model AD and give the corresponding manifold mor- 
phologies in Fig. [6l Some notable differences are clear (compare 
Fig. [6l of this paper with fig. 4 of Paper III), the most important 
one being that there are many more spirals here than for the previ- 
ously used rotation curve. This is in agreement with the results and 
discussions made so far, since higher Qt,L ± values lead to spiral 
morphologies and lower ones to an rRi morphology (Paper III and 
Sect l3.4l of the present paper). 

We also tried another axisymmetric component, such that the 
peak of the circular velocity curve is nearer to the galactic centre, 
i.e. more adequate for early type galaxies, and we denote the corre- 
sponding model by AE. This introduced relatively little change to 
the ring shapes, as will be discussed in the following sections. 



5.2.3 Outer ring axial ratios 

Analysing the results from his CSRG catalogue, lButaUl995h finds 
that the outer rings have a significantly oval intrinsic shape, with 
an average intrinsic axial ratio of < q Q >= 0.87 =b 0.14. For his 
SB sub-sample alone, he finds < q Q >= 0.82 ± 0.07. He also 
finds that the outer rings are rounder in SA than in SB galaxies, 
with an average intrinsic axial ratio of 0.97 ± 0.21. Applying the 
same procedure to the Southern G alaxy Catalogue (SGC) he finds 
< q Q >= 0.83±0.01 (lButall995h . while in a previous study with a 
somewhat different procedure he finds that the preferred axial ratio 
range is 0.7 - 1 for SB galaxies and 0.8-1 for SAs and SB As 
teutall 19861) . There is thus general agreement between the various 
estimates of the observed outer ring shapes. 

Amongst the roughly 300 models discussed in Sect. 3 of Pa- 
per III and in Sect. l5.2.2l above, 180 have an Ri or R[ morphology. 
This number is sufficiently large for us to be able to make some 
analysis. We therefore measured the major and minor axes of the 
outer ring for this sub- sample and plot their ratio as a function of 
Qt,L ± in Fig. [7] We have too few R2 models for any reasonable 
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Figure 7. Relation between the strength of the bar in the corotation re- 
gion, measured by Qt,Li> and the axial ratio of the outer R\, or R' x ring 
for our models. As in Paper III, different colours stand for different mor- 
phologies: green for R\ and orange for R' ± . The upper panel includes only 
morphologies for models A, D and BW and shows a strong correlation, in 
the sense that galaxies with stronger bars have also more elongated outer 
rings. Results from models A are given by filled diamonds, from models D 
by asterisks and from models BW by filled triangles. In the lower panel we 
add more models. Models AD are plotted with a filled circle and models 
AE with an open one. Models with lenses are plotted with an open square 
and models from Table|2]with an x sign. The horizontal lines give statistical 
information on the observed outer ring axial ratios, namely the mean (solid 
line) as well as the one a (dashed) and the two a (dotted) deviations from 
it, as calculated for barred galaxies from the CSRG (Buta 1995|). The gre y 
area shows the region of values found from the SCG survey (Buta 1986). 
Note that the outer ring axial ratio values we find from our models are in 
the right range. 



analysis and this morphology is also more rare in real barred galax- 
ies. Finally we can not compare the properties of R1R2 rings with 
observation, because the measurements are not uniquely defined 
sinc e it is not cle ar whether one should measure the Ri , or the R2 
part. iButal (1 19951) measures both and then takes an average of the 
values weighted according to the strength of each feature. We can 
not do this for our models because the manifold calculations we 
presented so far can give us information on the existence and shape 
of a feature, but not on how much material is trapped by it, the lat- 
ter depending on the evolutionary history of the galaxy. We thus 
have no information on how strong the Ri and R2 are. In order to 
obtain this information, we would need to follow the trapping by 
the manifold during the disc formation and evolution stage. A sim- 
plified version of such a calculation has been given in Sect. 14. 11 but 
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a more realistic version would necessitate considerable knowledge 
on how the disc and its main substructures were formed. 

The horizontal lines in Fig. [7] give the average value and the 
one an d two sigma deviations from it for barred galaxies from the 
CSRG (But all995h . T he grey area shows the region of values found 
from the SCG survey teutalll986h . We note that the values we find 
from our models are in the right ballpark. In fact, all values are 
within two sigma from the observed mean and most models are 
within one sigma from it. Adding results from other types of mod- 
els (with different axisymmetric components, or lenses, or stable 
Li and L2 Lagrangian points - see Sect. l5.2.2t does not change in 
any qualitative way these results (lower panel of Fig. [7]). Concern- 
ing the outer ring shape, there is thus very good agreement between 
theory and observations, which holds for all the very wide range of 
potentials we used. 

The upper panel of Fig. [7] shows a tight correlation, in the 
sense that galaxies with stronger bars have also more elongated 
outer rings. Most of the points on this plot correspond to A models, 
with some D and BW models. In fact, it would have been rather dif- 
ficult to assert the existence of this correlation from the D or from 
the BW models separately, because their Qt,L x values are rather 
clustered. Nevertheless, both the D ad BW models fall well on the 
correlation outlined by the A models, thus reinforcing it. 

Let us now check whether this correlation still holds for a yet 
wider set of model potentials (lower panel of Fig. [3, namely all 
those introduced in Sect. 15.2.21 All models with lenses fall neatly 
on the correlation. Most of models with a AD or AE axisymmet- 
ric component fall also on or very near the correlation. Only four 
points (out of 180) are considerably off. Furthermore, three of these 
correspond to marginal cases, i.e. cases which could be classified 
either as R[ , or as spiral. Thus we can conclude that, although the 
correlation is widened and could be called a trend, it is certainly 
still present, even with the very wide variety of the models used. 
This figure also gives an indication that, for this trend at least, the 
axisymmetric part of the potential, i.e. the circular velocity curve, 
makes more difference than the analytical expression used for the 
bar potential. Considerable more work is, however, necessary in 
order to assess this. 

On the other hand, models with stable L\ and L2 Lagrangian 
points definitely stand apart from the rest and are off the correlation 
although they have realistic values of the outer ring axial ratio. It is 
unclear at this point whether this argues only that this correlation 
will not be followed by models with stable L\ and L2 Lagrangian 
points, or, more strongly, that the correlation will be followed only 
by certain types of model potentials. 

It would be interesting to check whether the correlation 
we find in Fig. [7] is also valid for observed R\ and R[ outer 
rings. It is in agreement with the average values found by iButal 
dl995h . since these show that barred galaxies have more elon- 
gated outer rings than non-barred ones. Unfortunately, it is not 
possible to make any more accurate tests with only the CSRG 
data, because it is necessary to deproject each galaxy individ- 
ually before measuring the ring axial ratios. Furthermore, it is 
necessary to use an appropriate measure of the bar strength. 
The correlation in Fig. [7] uses Qt,L 1 which measures the non- 
axisymmetric forcing at corotation. Any such measure at this ra- 
dius, or somewhat beyond should be appropriate. Measures at 
smaller radii, however, may not be. To check this further, we re- 
plotted the data of the lower panel of Fig. [71 now using Qb as a 
bar strength measure. This quantity i s very often used in both theo- 
retical and observational studies ( e.g.lButa. Block & Kna pen 2003 ; 
iButa. Laurikainen &Sald 120041 : iLaurikain erTsalo & Butal 120041 : 



Figure 8. Outer ring axial ratio as a function of the bar strength. The latter is 
measured by Qb (see text for definition). Models A, D and BW are plotted 
by asterisks. The remaining symbols, and the colours, lines and shading are 
as in the lower panel of Fig. [7] 



lLaurikainen et allE004l : lDurbala et alll2009l and references therein). 
It relates to the maximum non- axisymmetric forcing and is defined 
as the maximum of Q t (r) in the bar region. 

The result is shown in Fig. [8] We have not included the results 
from the Dehnen models, because this ad hoc potential makes it 
difficult to measure Qb for the parameters we are using. Even with- 
out these points, however, there is no clear correlation. The best 
correlation is obtained if we look at models A alone and exclude 
all others. For this case, the rank correlation coefficient is -0.81, 
while the same models when plotted as a function of Qt,L x give 
a rank correlation coefficient of -0.94. Thus the correlation is less 
tight, but the trend itself is not lost. On the other hand, if we include 
also the BW models, then the rank correlation coefficient becomes 
-0.001 for Qb and -0.95 for Qt,L 1 , so that the correlation is quite 
strong for Qt,L 1 and is totally lost for Qb. Thus, whether it is possi- 
ble to use Qb instead of Qt,L x depends on the comparison at hand. 
In general, for properties related to the bar or to the inner manifolds 
Qb should be a good quantity to use. For the properties of the spi- 
rals and outer rings, however, the use of Qb may be hazardeous. In 
the absence of any further information it is safer and preferable to 
avoid the use of Qb for comparison with observations related to the 
spirals and rings. Nevertheless, Qt,L 1 is a more difficult quantity 
to obtain from observations than Qb, since it presumes an albeit 
rough knowledge of the Lagrangian radius tl- In Paper V we dis- 
cuss how this quantity can be obtained with the help of our theory 
and how reliable this estimate can be in various cases. One can thus 
conclude that such comparisons with observations should be done 
with great care since the value of Qt,L x may contain some error 
due to a wrong estimate of t\l, while Qb may give wrong estimate 
of Qt,Li and therefore of the quantities that depend on it. 



5.2.4 Inner ring axial ratios 

IButal dl995h analysed the inner ring data from the CSRG in a 
similar way as for the outer rings. He finds that the inner rings 
have a significantly oval intrinsic shape and in fact are on average 
more elongated than outer rings. For his SB sub- sample he finds 
< q >= 0.81 ± 0.06, where < q Q > is the average intrinsic axial 
ratio. As discussed above for outer rings, he can not seek correla- 
tions of the ring axial ratio with galaxy, or bar parameters because 
this would necessarily imply deprojecting each galaxy individually. 
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Figure 9. Relation between the strength of the bar near corotation, mea- 
sured by Qt,Lx » an d the axial ratio of the inner ring for our models. There 
is a strong correlation, in the sense that galaxies with stronger bars have 
also more elongated inner rings. The horizontal lines and the gray areas 
give some statistical information on the corresponding observed quantities 
(see text). Colours and symbols are as in Fig. [7] upper and lower panels, 
respectively. 



He does, however, find that inner rings are rounder in SA than in 
SB galaxies with an average intrinsic axial ratio of 0.9 2 ± 0.09 for 
the former. Applying the same procedure to the SGC. lButal (l995) 
finds < q Q >= 0.82 ±0. 02, while in a previous study with a some- 
what different procedure iButal dl986h finds that the axial ratios are 
uniformly spread in the range 0.60 - 0.95 for SB galaxies and 0.85 
- 1 for SAs. 

Fig. [9]gives the same information as Fig. [7] but now for inner 
rings. Results from model A (upper panel) show a very tight cor- 
relation between the inner ring axial ratio and the Qt,L ± measure 
of the bar strength. Models with lenses also lie on the same cor- 
relation. Models D and BW (upper panel) also show a correlation, 
which is followed by the models from Table[2] The regression lines 
of these two correlations, however, do not coincide, the former be- 
ing much steeper than the latter. 

A further point which is clear from the upper panel of Fig. [9] 
is that some of the models have an inner ring shape which is much 
more elongated than the observations. Several arguments need to be 
taken into account here. It should be kept in mind that all the biases, 
both observational and theoretical, tend, as discussed in Sect. l5.2.fl 
to increase this discrepancy and their cumulated effect could be 
considerable. Amongst the most important ones is the fact that, as 
we already mentioned, the inner manifold branches can correspond 



either to the ring or to the outer parts of the bar. For models A, AD 
and AE where the bar outline is known, it is possible to distinguish 
between the two. Unfortunately this is impossible for models D and 
B W and not straightforward for models from Table[T]or Tabled We 
examined all our A, AD and AE models and eliminated those for 
which the inner manifolds are more likely to correspond to the bar 
than to an inner ring. This exercise is not trivial, because our cal- 
culations are not self-consistent, i.e. the gravity from the manifold 
itself is not taken into account. If the manifold locus is somewhat 
outside, but near the bar outline, then in real galaxies, or in self- 
consistent calculations, the bar outline could somewhat change due 
to the gravitational pull from the manifold and, provided the dis- 
tance between the two is small, the two could merge. This is some- 
thing we can not take into account with our calculations. We thus 
eliminated only the cases where we estimated from our non- self- 
consistent simulations that the manifold pertained to the bar more 
than to the inner ring. We then plot the remaining models, including 
all possible families, in the lower panel of Fig. [9] 

Note that most of the models are in the area outlined by the 
observations, but a handful are still clearly outside. We examined 
by eye these cases and found that they are borderline, i.e. perhaps 
a self-consistent calculation would lead us to classify them as con- 
tributing to the outer part of the bar, rather than to the ring. 

We also examined which potentials increase this discrepancy, 
i.e. make more elongated inner rings, and which alleviate it. First 
considering the rotation curves, we see that AD rotation curves 
create more elongated rings, i.e. increase the discrepancy between 
models and observations. This can be understood since in AD mod- 
els the axisymmetric force is relatively smaller over some of the 
relevant radial extent than in A models and this, according to our 
previous results, should lead to more elongated rings, as is indeed 
found. For models AE the difference is less clear-cut. We can, how- 
ever, see that for the largest Qt,L t values, the AE have less elon- 
gated rings than the A models, thus alleviating the discrepancy. 
This can be understood by the fact that, for the most elongated 
rings, the inner manifold branches reach radial regions where the 
axisymmetric force is relatively smaller in the AE models than in 
the A ones. This would then lead to less elongated inner rings in 
the cases with largest Qt,L 1 , as indeed found in our calculations. 

Th e problem of too elongated inner rings was already encoun- 
tered bv lSchwarzUl985h , who attributed it to the potential he was 
using. To solve it he proposed to include a lens component. Fig. [9] 
and Table Q] show that this is indeed a possible solution. For exam- 
ple, the inner ring axial ratio of one of the A models, which without 
a lens is equal to 0.51 (i.e. nearly three sigma below the average), 
becomes 0.75 (i.e. within one sigma from the mean) if half of the 
bar mass is turned into a lens. This of course can not be the only 
solution, since not all galaxies with inner rings have lenses. 

The above argue convincingly that the axial ratio of the in- 
ner ring depends on the potential used, stronger relative non- 
axisymmetrc forcings leading to more elongated rings. The results 
we get with our potentials are in relatively good agreement with 
observations and an excellent agreement could be reached if the 
non-axisymmetric forcings in the models were somewhat lower in 
the region covered by the inner ring. 

5.2.5 Statistics on the major axis ratios 

Analysis of the ratio of outer to inner ring sizes in observed galaxies 
has shown that the c orresponding histogram is rather concentrate d 
to values around 2 (I Athanassoula et al |[l982l : lButalll98fl Il995h . 
This is a clear-cut result and it would be nice to test our models 
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Figure 10. Ratio of the major axis of the outer ring to the major axis of 
the inner ring as a function of Qt,L 1 for all our models. Green symbols 
denote measurements from R\ models, and orange ones from R f ± models. 
The symbol used depends on the model potential. Results from models A 
are given by diamonds, from models D by asterisks and from models BW 
by filled triangles. Models AD are plotted with a filled circle and models 
AE with an open one. Models with lenses are plotted with an open square 
and models from Table|2]with an x sign. 
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Figure 11. The effect of a change of the model circular velocity curves on 
the ratio of the outer to inner major axes. The bar parameters do not change, 
so that each point shows the difference between two models with the same 
bar, but different circular velocity curves. On the abscissa we plot the dif- 
ference between Qt,L ± for model AD and the same quantity for model A. 
On the ordinate we plot the difference between the ratio of the outer to inner 
major axes, D /Di , for models AD and the same quantity for model A. We 
note that the change from A to AD increases both the value of Qt,L 1 an d 
of D /Di, thus bringing the latter measurements nearer to the observed 
values. 

with it. However, as already discussed in detail in Sect. 15.2.11 this 
is a very difficult comparison to make, since the ratio in question 
is strongly affected by errors and biases due to both difficulties in 
measuring the observed and the model ring sizes and to the fact that 
measurements in observations have not been done in individually 
deprojected galaxies. Finally, it should be remembered that we can 
compare the range of values covered by observations, but we can 
not compare the corresponding distributions, since the models we 
have presumably do not cover the parameter space in the same way 
as the observed galaxies. In fact the observed distribution of some 
of the relevant parameters is still largely un known, althou gh consid- 
erable progress has been made for others dGadottill2009h . We will, 



nevertheless, attempt a general discussion and comparison with ob- 
servations and try and outline which types of potentials give the 
best (or the worst) results. 

Fig. [TOl shows the ratio of the outer to inner ring major axes, 
D /Di, plotted as a function of Qt,L 1 - It is clear that although a 
large number of our models have values in the same ballpark as the 
observations, many have considerably smaller values. Fig. [Tol also 
shows a trend which, although it is not sufficiently tight to be called 
a correlation, is still very clear. Namely, the ratio D Q / Di is smaller 
in cases with smaller Qt,L 1 values. Fig.QIJalso gives valuable in- 
formation on which types of potentials are in better agreement with 
the observations. 

All models with an R[ morphology have a ring ratio in agree- 
ment with observations, independent of the type of potential used. 
The same is true for nearly all models of the AD type, as well as 
for models whose L± and L2 Lagrangian points are stable. All these 
give results which are roughly in the right range of values. On the 
other hand, many of the remaining cases give smaller values of the 
axial ratio, i.e. outside the desired region. Such results come from 
models with an rR± morphology with potentials A, or AE. 

To make the effect of the circular velocity curve clearer we 
compare in Fig. Qj] models A to models AD. For this we consider 
pairs of models having the same bar parameters, but one having an 
A circular velocity curve and the other an AE one, and calculate 
the difference in Qt,L ± and in the ratio of outer to inner ring major 
axes, D /Di. We see that both quantities are larger for AD models 
and that the difference can be considerable. Thus such a change 
would bring our D Q / Di in better agreement with observations. 

This clearly underlines the effect of the properties of the 
model potentials on quantitative measures and fully justifies the 
corresponding arguments advanced in Sect. 15.2.11 It should not, 
however, be interpreted as arguing that rotation curves should be 
slightly decreasing at and beyond corotation. It argues instead for 
somewhat larger Qt,L ± values in that region, something that could 
be also due to a somewhat larger non-axisymmetric forcing. This 
would be the case in self-consistent treatments, since the potential 
from the non-axisymmetric rings would be included, and also if the 
bar model potential dropped less abruptly with radius that the n = 1 
Ferrers models we use. Thus, such an increase seems very natural 
and it would easily lead to values of the ring size ratios in good 
agreement with observations. 



6 SUMMARY 

In previous papers (Papers I, II and III) we proposed a theory to 
explain the formation and properties of rings and spirals in barred 
galaxies. In this paper we compared their morphological properties 
to that of observed ones. Our theory can produce both spirals and 
inner and outer rings. Furthermore we get all three observed types 
of outer rings , i.e. R±, R2 and R1R2 rings and pseudo-rings, and 
no other morphologies. This means that all outer rings produced 
by our theory have a qualitatively realistic morphology and that 
all observed morphologies are qualitatively reproduced. We also 
reproduce the dimples often seen clearly in Ri rings. Furthermore, 
we predict that if the non-axisymmetric forcing is relatively low the 
resulting morphology will be an R± ring, while if it is stronger, it 
will be a spiral. 

Spirals due to this theory start off from a location not far from 
the end of the bar and are trailing, in good agreement with obser- 
vations. They are in the vast majority of cases two armed, although 
it is possible, if specific potentials are used, to have arms of higher 
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multiplicities. Their shape is also realistic and we made a success- 
ful quantitative comparison with the arms of NGC 1365. They are 
able to reproduce the "fall back" of the arm towards the more in- 
ner region, seen in NGC 1365 and in many other barred spirals, 
something which density wave theory could not achieve. We also 
predict that stronger non-axisymmetric forces in the region around 
and somewhat beyond corotation will result in more open spirals. 

Quantitative comparisons are less straightforward to make, not 
because they imply tighter constraints, which they do not, but be- 
cause they involve not only the theory, but also the potential to 
which it is applied. We are confident that our theory is dynamically 
correct, and the comparisons with observations made here and in 
Paper V make us optimistic about its applicability to galaxies with 
spirals or rings. On the other hand, we are aware of the shortcom- 
ings of our potentials. First and foremost, our analysis is not self- 
consistent, so that we neglect the self-potential of any structure they 
create. Furthermore, the Ferrers' potentials, although the best of all 
the available analytical bar potentials, have a number of shortcom- 
ings, the most important one for our particular application being 
that the force they create diminishes too rapidly at large radii. The 
other potentials we used, although they lack this particular draw- 
back, are ad hoc, and therefore not necessarily realistic. The above 
remarks should be kept in mind when assessing our quantitative 
comparisons. Our aim in such cases was to find whether the poten- 
tials we used give results in agreement with observations, or, in case 
of disagreement, whether any obvious and straightforward changes 
of the potential can bring agreement. 

For the outer ring our analysis concentrates on R± cases, for 
which we have sufficient models to reach safe conclusions. With 
all the potentials used so far, the outer ring shape is in good agree- 
ment with observations. Most models are within one a from the 
observed mean and all are within two a. We tried different changes 
of the potential - such as different rotation curves, adding a lens, or 
stabilising the L± and L2 Lagrangian points (see Paper III) - and 
this result remained unchanged. There is thus excellent agreement 
between model and observed outer ring shapes. 

For the inner ring, we found several models with inner ring 
shapes more elongated than observed. Examining by eye these 
cases, we found that for most, if not all, the inner manifolds cor- 
respond more to the outer regions of the bar and not to the inner 
ring. We can thus conclude that there is a fair agreement between 
theory and observations regarding the inner ring axial ratio. 

Finally, the ratio of outer to inner ring major axes gives good 
agreement for some types of models, but not for all. Thus all 
pseudo-ring cases (rR[), independent of the model potential, give 
good agreement. This is the case also for both rR± and rR[ mor- 
phologies in models with a circular velocity curve showing a slight 
decrease at or somewhat beyond corotation, for models with stable 
Li and Li Lagrangian points and for a number of other models. On 
the other hand, agreement for cases with an rR\ morphology and A 
or AE rotation curves (with or without a lens) often have too small 
a value of this ratio. It is thus clear that the underlying potential 
plays an important role. Our calculations suggest that models such 
that Qt,L 1 is somewhat decreased within corotation and increased 
outside it would give good agreement for all ring size measure- 
ments, namely the inner and outer ring shapes (ratio of minor to 
major axes) and the ratio of the outer to inner ring major axes. It is, 
however, beyond the scope of this paper to find a model achieving 
this at all relevant radii, thus includes the best parts of all the mod- 
els considered, particularly since this exercise may prove futile if 
self-consistency is taken into account. 

We find a trend between the ring axial ratio and the relative 



amplitude of the non-axisymmetric forcing in the region around, 
or immediately beyond corotation, in the sense that stronger bars 
are linked with more elongated rings. This is true for both inner 
and outer rings. There is also a trend between the ratio of outer to 
inner ring major axes and this measure of non-axisymmetric forc- 
ing. If we focus on only one model family, or a couple of them, 
these trends can become clear correlations, but including all mod- 
els broadens them, sometimes substantially. 

In Paper V, the next and last paper of this series, we will 
present a number of other comparisons with observations and pre- 
dictions, including kinematics, bar shape, radial drift and abun- 
dance gradients. We will discuss pattern speed prediction and give 
a global appreciation of the application of our theory to observed 
spirals and rings. We will also consider the effects of time evolution 
and of self-consistency. Finally we will discuss other theories and 
simulations and how they link to this work. 
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